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(SI . ABSTRACT 

We have characterized the physical properties (electron temperature, density, 
. metallicity) of the ionized gas and the ionizing population (age, metallicity, pres- 

ence of WR stars) in the Lynx arc, a HII galaxy at z =3.36. The UV doublets (CHI], 
Silll] and NIV) imply the existence of a density gradient in this object, with a high 
density region (0.1-1.0 x 10 5 cm~ 3 ) and a lower density region (<3200 cm -3 ). The 
temperature sensitive ratio [OIII]AA1661,1666/A5007 implies an electron temperature 
. T e =173Q0^ 70 g K, in agreement within the errors with photoionization model predic- 

(n^| ' tions. Nebular abundance determination using standard techniques and the results 

\ from photoionization models imply a nebular metallicity of O/H~10±3% (O/H)0, 

C*~) ■ in good agreement with Fosbury et al (2003). Both methods suggest that nitrogen is 

0^ ' overabundant relative to other elements, with [N/0]~2. 0-3.0 x [N/O]©. We do not 

, find evidence for Si overabundance, as Fosbury et al. (2003). 

>J ■ Photoionization models imply that the ionizing stellar population in the Lynx arc 

has an age of Ss5 Myr. If He + is ionized by WR stars, then the ionizing stars in the 
. Lynx arc have metallicities Z star >5% Zq and ages ^2.8-3.4 Myr (depending on Z star ), 

when WR stars appear and are responsible for the He +2 emission. However, alternative 
1~ \ excitation mechanisms for this species are not discarded. Since the emission lines trace 

the properties of the present burst only, nothing can be said about the possible presence 
' of an underlying old stellar population. 

, The Lynx arc is a low metallicity HII galaxy that is undergoing a burst of star 

■ formation of ^5 Myr age. One possible scenario that explains the emission line spec- 

, £^ ' trum of the Lynx arc, the large strength of the nitrogen lines and the He +2 emission 

, is that the object has experienced a merger event that has triggered a burst of star 

5_i ■ formation. WR stars have formed that contribute to a fast enrichment of the ISM. 

As Fosbury et al. (2003), we find a factor of ^10 discrepancy between the mass 
of the instantaneous burst required to power the luminosity of the H/3 line and the 
mass implied by the continuum level measured for the Lynx arc. We discuss several 
possible solutions to this problem. The most likely explanation is that gas and stars 
have different spatial distributions so that the emission lines and the stellar continuum 
suffer different gravitational amplification by the intervening cluster. 

Key words: cosmology: observations - galaxies: abundances - galaxies: high redshift 
HII regions - stars: formation 



1 INTRODUCTION 

H II galaxies are dwarf emission-line galaxies undergoing a 
burst of star formation. They are characterized by strong 
and narrow emission lines originated in a giant star-forming 
region which dominate their observable properties at opti- 
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cal wavelengths (e.g. (Terlevich et al. 2004)). Most are blue 
compact galaxies (BCGs) . They have very low metallicities, 
high rates of star formation and a very young stellar con- 
tent. Many are compact and isolated. One of the reasons 
why these objects have attracted significant attention is the 
possibility that they are very young galaxies in the process 
of formation. This possibility, however, has been challenged 
since evidence for an old (several Gyr) stellar population has 
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been found in numerous BCGs (see (Kunth & Ostlin 2000) 
for a review). Therefore, a model with a single, instanta- 
neous burst of star formation does not seem appropriate to 
describe these galaxies and a succession of short starbursts 
separated by quiescent periods seems more likely ((Terlevich 
et al. 2004)). 

Even in the case if HII galaxies are not primeval, little 
chemical evolution has happened and they provide impor- 
tant information about how galaxies form and evolve, as well 
as about the process of star formation in low metallicity en- 
vironments. 

The most distant HII galaxy in the catalogue compiled 
by Terlevich et al. (1991) is at redshift z ~0.31, while the 
vast majority are at z <0.1. The proximity of these galaxies 
allows studies of their structure, metal content and stellar 
population with high sensitivity and precision, difficult to 
achieve for high redshift star forming galaxies. Still, finding 
high redshift HII galaxies is of great interest to investigate 
their abundance patterns and stellar properties at earlier 
epochs, compare them with the nearby counterparts and ob- 
tain further information about the formation and evolution 
of galaxies using objects that have undergone little chemical 
evolution. 

The Lynx arc is a star forming object at z=3.36. It was 
discovered during spectroscopic follow-up of the cluster RXJ 
0848+4456 (z=0.57) from the ROSAT Deep Cluster Survey 
((Holden et al. 2001)). The arc shows a very red R-K color 
and very strong optical and UV (rest frame) narrow emission 
lines. From the analysis of HST WFPC2 images and Keck 
optical and near infrared spectroscopy, Fosbury et al. (2003, 
FOSB03 hereafter) concluded that the arc is an HII galaxy 
magnified by a factor of ~ 10 by a complex intervening clus- 
ter environment. The authors concluded that the continuum 
is mostly nebular. 

By means of photoionization modeling, (FOSB03) 
showed that the spectroscopic properties of the Lynx arc are 
all consistent with a simple HII region model of a cluster of 
~10 massive stars, characterized by a very high T e // ~80 
000 K, a high ionizing parameter (U~0.1) and low nebu- 
lar metallicity Z~ 5% Z . Indirect arguments suggest that 
the ionizing stellar population could have much lower abun- 
dances. Using simple photoionizaton models, the authors 
conclude that the spectroscopic properties of the Lynx arc 
are consistent with those of a metal-poor nebula ionized by a 
cluster of primordial (i.e. metal-free) stars. In this scenario, 
the overabundance of Si implied by the models is explained 
as due to enrichment by past pair-instability supernovae, 
requiring stars more massive than 120 M . 

We investigate in this paper the possibility that normal 
(i.e. non primordial) stars are responsible for the excitation 
of the gas in the Lynx arc. Our first goal is to constrain the 
physical properties of the gas: electron temperature, density, 
metallicity. We will first use standard spectroscopic tech- 
niques of nebular analysis (e.g. (Aller 1984)). These results 
will be compared with those from detailed photoionization 
modeling and those obtained by (FOSB03). 

The second goal is to characterize the ionizing stellar 
population: age, metallicity, stellar mass of the burst, possi- 
ble presence of WR stars. This will be achieved by means of 
detailed photoionization models, whose objective is to find 
an ionizing stellar population able to explain all the spec- 
troscopic properties of the Lynx arc (line ratios, equivalent 



widths, intensity of the continuum). In addition, such mod- 
els set tight constrains on the nebular properties, which must 
be consistent with the results from the standard techniques 
mentioned above. 

We present in §2 the results of the nebular analysis and 
the photoionization models. These results will be discussed 
in §3. Summary and conclusions are presented in §4. 



2 ANALYSIS 

2.1 Nebular properties 

In this section we will set constrains on the electron tem- 
perature, density and the nebular abundances of the ionized 
gas in the Lynx arc. 

All the analysis and discussion presented below are 
based on spectra obtained with different slit widths and 
orientations (see (FOSB03) for a detailed description of 
the data set and reduction techniques). For this reason, 
(FOSB03) estimated carefully scaling factors that were ap- 
plied in order to derive a consistent calibration. We will 
therefore assume that the errors introduced by calibration 
uncertainties are negligible. 

The results are also dependent on whether reddening is 
present or not. We have ignored this effect. The very strong 
UV lines and the lack of evidence for dust reddening from 
the photoionization models and the fit to the continuum 
shape ((FOSB03)) makes as confident that such assumption 
is reasonable and our results are not seriously affected by 
it. Our own work presented in this paper will show that the 
line ratios are consistent with no dust reddening. 



2.1.1 Electron temperature 

The [OIII]4363/[OIII]5007 ratio has been traditionally used 
to measure electron temperatures in ionized nebulae. Unfor- 
tunately, the [OIII]A4363 line is outside the observed spec- 
tral range. We have used the UV [OIII] doublet instead, since 
the [OIII]AA1661,1666/[OIII]A5007 ratio is mostly sensitive 
to electron temperature. 

We show in Fig.l 

the dependence of the [OIII]AA1661,1666/[OIII] A5007 ratio 
with electron temperature The position of the Lynx arc (the 
value of the ratio is 0.074±0.005, see Table 1) is also shown 
as a filled circle. It implies T e ~17300±^ K - This value is 
somewhat higher than that predicted by the photoioniza- 
tion models presented in §2.2 which give T e =16200±500 K, 
although taking errors into account, the discrepancy is not 
important. This further supports that dust reddening effects 
are negligible, otherwise, the temperature derived from Fig. 
1 would be a lower limit and large discrepancies with the 
photoionization model predictions would arise. This tem- 
perature is much lower than that derived by (FOSB03) and 
(Binctte et al. (2003)) by means of photoionization model- 
ing (~20000 K). The discrepancy is due to a combination of 
the very hard ionizing continuum (very energetic electrons 
are therefore released in the ionization processes) and high 
U (ionization parameter) value used by these authors. 
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Figure 1. Dependence of the [OIII] AA1661,1666/[OIII] A5007 ra- 
tio with electron temperature T e .The position of the Lynx arc is 
indicated with a solid circle. This diagram implies T e ~17300 4 ^?q 
K for this object. 
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are also emitted within this region. There is a zone of lower 
ionization level, where the C +2 , Si +2 and N +2 lines are emit- 
ted and which has n e <3200 cm -3 . 

Usually, the optical estimates of electron densities of 
H II regions, based on forbidden emission line ratios, are 
in the range 100 to 1500 cm -3 . However, radio techniques 
have revealed the existence in the Galaxy of very small and 
dense (ri4 >10 4 -10 5 cm -3 ) objects, the so-called compact or 
ultracompact HII regions (e.g. (Habing & Israel 1979)). 

Density gradients in several giant and galactic HII re- 
gions have been observed (e.g. (Copetti et al. 2000)). In 
most of these cases, the spatial variation of density may be 
interpreted as a radial gradient with the density decreasing 
from the centre to the edges. A good example is the Orion 
nebula. Based on measurements of the ratio [O II]3729/3726, 
Osterbrock & Flather (1959) showed the existence of a steep 
radial density gradient in the Orion Nebula, with the elec- 
tron density decreasing from 1.8xl0 4 cm -3 in the centre to 
2.6xl0 2 cm -3 near the edge. 



2.1.3 Nebular abundances 

Photoionization models by (FOSB03) imply that the nebu- 
lar abundances in the Lynx arc are ~5% solar within a fac- 
tor of 2. We set here tighter constraints by means of using 
standard techniques of nebular abundance determination. 

The following expressions were used to calculate the 
ionic abundances of several elements ( Aller 1984) , assuming 
the low density regime. These expressions will then be used 
to calculate total element abundances. 



Table 1. Electron densities measured for the Lynx arc using the 
CIII], Silll] and NIV] UV doublets. 



2.1.2 Electron density 

There are three UV line doublets we can use to estimate 
the electron density in the Lynx arc: CIII] AA1907,1909, 
SiIII]AA1883,1892 and NIVAA 1483, 1487 (see (Keenan, 
Feibleman & Berrington 1992), (Keenan et al. 1995)). The 
results shown in Table 1 correspond to T e ~ 10000-20000 K 
range, where the ratios are quite insensitive to temperature. 
The upper limits take the errors into account (i.e. the densi- 
ties shown correspond to the minimum possible value of the 
line ratio). 

The Silll] doublet density determination has to be con- 
sidered carefully. We have noticed that the predictions by 
Keenan, Feibleman & Berrington (1992, see Fig. 2) are ac- 
tually inconsistent with the results of our photoionization 
models (§2.3). As an example, for some of the n e = 1000 
cm" 3 models presented in §2.3, SiIII]1883/1892=1.8, which 
is above the maximum possible value (~1.5) predicted by 
those authors for this ratio. For models with n e =10 cm -3 , 
we find that the ratio reaches a value of 2 in some cases. 
We also find that the ratio varies for a fixed density (be- 
tween 1.45 and 1.8 for n e =1000 cm -3 ), depending on other 
parameters such as the spectral energy distribution (SED) 
properties. 

The range of densities implied by the different doublets 
suggests that there is a gradient in the Lynx arc. The highly 
ionized region responsible for the emission of the N +3 lines 
has n e ~[0.1-1.0]xl0 5 cm -3 . The C+ 3 lines (1549,1551 A) 
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For silicon we have used the expression: 
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Si +2 nl oo f o.2 om/t SUII] 1883, 1892 
C+ 2 6 CJJJ] 1906, 1909 

(from (Garnett et al. 1995)) 

Garnett et al. (1995) showed that this expression is a 
good approximations within the temperature range 10000- 
20000 K and for the n e values expected for the Lynx arc, 
well below the critical densities of all the lines considered 
above. 

The total element abundances are then calculated as 
follows: 

O 0+ 2 Ne Ne+ 2 

H ~ ~H+ IT ~ H+ 

In both cases, we assume that most oxygen and neon 
atoms are twice ionized and therefore the ion abundance ra- 
tio is a good representation of the element abundance ratio. 
This assumption is supported by the high level of ionization 
of the gas implied by the high [OIII]/H/3 and [OIII]/[OII] 
ratios of the Lynx arc. It is also confirmed by our photoion- 
ization models (see §2.2) which show that ~99% of O and 
Ne are in the +2 and Ne +2 forms. 

For the other elements we use: 

c c+ 3 + C+ 2 

H H+ 

N N+ 4 + N+ 3 + N +2 N +3 + N+ 2 
H H+ H+ 

(as Table 1 shows, the abundance of N +4 in negligible 
compared with that of N +a and N +2 ) 

We assume that the low ionization species N + , N°, C + 
and C° are not present, which is expected due to the high 
level of ionization of the gas. It is also confirmed by our 
photoionization models (see §2.2), that also predict negligi- 
ble abundance for ionization species higher than the ones 
considered above. 

To calculate the abundance of Si we need to account 
for the presence of different ions of this element using an 
ionization correction factor ICF: 

si_s^ x MSi+ 2 ) t Si+ 2 

c - c+2 x 1 X (C+ 2 ) 1 ~ C+ 2 1 ^ rsi ' c 

Garnet et al. (1995) predicted the values of ICFsi/c 
for different X(0 + ) (see Fig. 1 in their paper). The scat- 
ter at large X(O +2 )(~0.99) values implies that ICF Si / C is 
in the range 2.5 - 5. Our photoionization models (§2.2) im- 
ply ICF Si /c in the range 3.2-3.4 (best models with no WR 
stars) and 2.3-2.9 (best models with WR stars). We there- 
fore assume lCFgi/c in the range 3.4 - 2.3 (accounting for 
different stellar metallicities/ages). 

The results of the calculations are shown in Table 2. Er- 
rors were calculated using the uncertainties on the line mea- 
surements presented in Table 3^. The errors on the abun- 
dances relative to the Sun are less than few per cent. C and O 

t Errors on line measurements provided by A. Humphrey, private 
communication. We have not added an error for CIV/H/3, since 
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Table 2. Abundances of some ions and elements in the Lynx arc 
for two different t values 1.6 and 1.7. The range given for Si cor- 
responds to the extreme values of the ionization correction factor 
ICF implied by our photoionization models (see text). The sec- 
ond value corresponds to models with WR stars. The two values 
given for oxygen correspond to the calculations performed using 
the [OIII]AA4959,5007 and the [OIII] AA1661,1666 doublets. The 
errors on the abundances relative to the Sun are less than few per 
cent. This only considers the errors on the ionic abundances and 
do not account for uncertainties on the Solar values. 

solar abundance values come from photospheric abundances 
of Allende-Prieto, Lambert & Asplun (2002, 2001), while N, 
Ne, Mg, Si and Fe are from (Holweger (2001)). Uncertainties 
on these values were not considered in the error calculation. 

Calculations were done for two temperature values 
T e =16 000 and 17 000 K (t=1.6 and 1.7) to account for the 
most likely T e range implied by the photoionization models 
(§2.2) and the [OIII] lines (§2.1.1). It is important to note 
that the electron temperature is different in high and low 
ionization zones of HII regions. To perform a more accurate 
abundance determination, this temperature variation should 
be taken into account. Most of the gas is ionized in the Lynx 
arc (~99 % of oxygen in the +2 form) and there is not a 
low ionization zone (i.e. a region where species such as [OI], 
[Oil], [SII], etc exist). We expect therefore that the electron 
temperatures measured using the [OIII] lines can be safely 
assumed to be representative of a large fraction of the neb- 
ula. Somewhat higher temperatures might be present in the 
most highly ionized region. In such case, the abundances 
presented in Table 2, are upper limits t. 

We conclude that the nebular metallicity in the Lynx 

this value is determined by the uncertainties on the reconstruction 
model of the line profile, rather than the line flux measurements 
((FOSB03)) 

t The small discrepancy between the O/H values calculated with 
the [OIII]AA5007,4959 and OIII]AA1661,1666 values (see Table 
2) could be due the fact that the excitation of the UV doublet 
requires somewhat higher electron temperatures than the optical 
doublet. 
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arc (defined as 0/H) is in the range ~10 - <20% Z . This 
is somewhat higher than the value derived by (FOSB03) of 
-5% Z©. 

In addition, we conclude that nitrogen is overabundant 
relative to oxygen by a factor of ~2. 0-3.0 compared with the 
solar abundance. 

Departures from the constant density case may produce 
variations in the derived ICF values and therefore, the el- 
ement abundances (Luridiana & Cervino 2003, Luridiana 
et al. 2004, in prep.). This does not affect oxygen seriously. 
However, the derived abundance of nitrogen can be off by up 
to a factor of 1.75 for their models with the steepest power 
law density gradient (index >1.5). It is therefore possible, 
that N is not overabundant and the large derived N/O is a 
consequence of assuming constant density. 

(FOSB03) concluded that silicon has to be well over- 
abundant relative to other elements, to be able to explain 
the strong Silll] lines detected in the Lynx arc. The reason 
for this discrepancy will be discussed in §2.2. 

Photoionization models in §2.2 will set further con- 
strains on the nebular abundances of all elements. 

2.2 Age and metallicity of the ionizing stellar 
population 

We have characterized the age and metallicity of the ioniz- 
ing stellar population in the Lynx arc by means of produc- 
ing photoionization models that best reproduce the emission 
line properties. Further constraints on the nebular physical 
properties and chemical abundances will also be set. 

The most outstanding spectroscopic properties of the 
Lynx arc that models should reproduce are (see (FOSB03)): 

• The large strength of the UV collisionally excited lines 
(NIVA1486, CIVA1550, OIIIJA1663, NIII]A1750, etc) relative 
to recombination lines such as H/3 

• The large [OIII]/H/?(=7.5) and the weakness of 
[OIIJA3727 relative to other lines ([OIII]/[OII] >30, 
[OII]/H/3<0.25 and [Neffl]A3869/[OII]>2.8)). I.e., the high 
ionization level of the gas. 

• The strength of the He +2 emission 
(HeIIA1640/H/3=0.11) 

• The extremely weak UV and optical continuum 

All the SEDs used in this paper were computed by 
Cervino, Mas-Hesse & Kunth (2004, hereafter CMK04) 
available at http://www.laeff.esa.es/users/mcs/SED. 
The isocrones were computed following the prescriptions 
quoted in (Cervino, Gomez-Flechoso & Castander (2001)) 
and the models assume a power law Initial Mass Function 
(IMF) with a Salpeter slope, a=2.35 ((Salpeter 1955)) in 
the mass range 2-120 Mq. We have adopted the evolution- 
ary tracks with standard mass loss rates by (Schaller et 
al. (1992); Schaerer et al. (1993a); Schaerer et al. (1993b); 
Charbonnel et al. (1993)) and the following atmosphere 
models: (Schaerer & de Koter (1997)) CoStar for main- 
sequence hot stars more massive than 20 Mq, (Schmutz, 
Leitherer & Gruenwald (1992)) for WR stars and (Kurucz 
(1991)) for the remaining stars. An instantaneous star for- 
mation law was assumed. We have also considered the in- 
fluence of the X-ray emission of the starburst in the SED 
by varying the efficiency of conversion of kinetic energy into 
X-rays. However, we have found that this effect does not 



produce significant differences in the predicted emission line 
spectrum at the young ages considered here. Therefore, we 
have ignored it and we present only the results for SED 
models without X-ray emission. 

We used the numerical code Cloudy ((Ferland et al. 
1998)) to predict the line ratios in different circumstances. 

Given the complete lack of information on the geometry 
of the gas and the relative distribution between stars and gas 
in the arc, we have assumed the standard spherical geometry 
often used in photoionization models of other extragalactic 
HII regions and HII galaxies. The distance between the ion- 
izing source and the gas ro was assumed to be 10 pc and we 
used constant density n=1000 cm -3 . The filling factor // 
was varied until finding the model that best reproduces the 
Lynx arc line ratios. Conclusions on the size and geometry 
of the ionized region cannot be extracted from our modeling 
since different values of //, n and ro might produce the same 
output spectrum. In addition, the assumption of spherical 
geometry is likely to be too simple, since the ionizing stars 
might be distributed in several clusters. 

We adopt log Q(H) = 55.20, as calculated by (FOSB03) 
from the H/3 flux and corrected for magnification. Q(H) is 
the ionizing luminosity in erg s _1 . Since the authors assumed 
covering factor CC=1 in these calculations, we have also 
adopted this value. 

We show in Tables 3 to 6 the results of the models 
that best reproduce the Lynx spectrum, together with the 
Lynx measurements. Although He +2 is a problem for some 
of these models, they are also shown because of the very 
good agreement with the rest of the line ratios. The He +2 
problem will be discussed in more detail in §3.1. We find 
that the nebular metallicity Z ne & that produces the best fit 
to the Lynx spectrum is ~10% Zq. This is in very good 
agreement with the results obtained in §2.1.3. 

Z„ e b=5% Zq models are rejected. The models that can 
produce [OIII]5007/H/3 J>7, require very high filling factor 
// £0.32 (i.e. high ionization level). The predicted electron 
temperature is ~20000 K, much higher than the value im- 
plied by the [OIII] lines (§2.1.1). In addition, due to the high 
T e and ionization level of the nebula, the high ionization UV 
lines such as CIVA1550 and [OIII1A1665 become too strong, 
with CIV/H/3 > 7. 

The opposite effect is found for models with Z ne b >15% 
Zq , which can also be rejected with confidence. Due to the 
higher metallicity, lower filling factors (i.e. ionization level) 
are required to produce the same [OIII] A5007/H/3. For those 
models with [OIII]5007/H/3 £7-8, the predicted T e is ~14 
000K, too low compared with the expected value (§2.1.1). 
In addition, because of the low ionization level and the low 
electron temperature, the high ionization UV collisionally 
excited lines become too weak, with CIV/H/3 

Therefore, all the models presented here assume 
Z ne 6 = 10% Z for all elements (except nitrogen, see below). 
The most discrepant line ratios will always be shown in bold 
characters. 

Zstars — 5% Z models 

The results of these models are shown in Table 3. The 
2 Myr(a) model assumes that all elements have an abun- 
dance of 10% the solar value. Notice that the NIII] and NIV 
lines are severely underestimated. This is the case for all 
models considered below, as long as N abundance is not en- 
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hanced. The discrepancy disappears if we assume that this 
element is overabundant by a factor of ~3. This result was 
also suggested by the nebular abundance determination in 
§2.1.3. Hereafter, this overabundance will be assumed in all 
models. 

Models with ages <4 Myr reproduce very well the Lynx 
line properties, except the He +2 emission. Ages of 5 Myr and 
specially older present strong discrepancies with the mea- 
sured line ratios, due to the drastic reduction of the ionizing 
photons as stars get older. We only show 5 and a 6 Myr 
models in Table 3 for illustration. 

The He +2 lines are always predicted too faint. The rea- 
son is that according to the stellar evolution models no WR 
stars are formed in such low metallicity environment (but 
see §3.1). 

Zstars = 20% Z Q models 

Models with ages ^3 Myr are in very good agreement 
with the Lynx spectrum (Table 4) , except for the He +2 lines 
which are predicted too faint. At ~3.4 Myr Wolf Rayet stars 
make an important contribution to the the hard ionizing 
radiation and the He +2 problem is solved. NV should be 
detected according to this model. However, given the uncer- 
tainties on the estimation of the upper limit for this line, 
the predicted ratio is likely to be well within the errors 
((FOSB03)) of the measured value. Ages of 4 Myr and older 
produce strong discrepancies with the measured spectrum. 
We show 4 and 5 Myr models in Table 4 for illustration. 

Tistars = 40% Z models 

Models with ages <3 Myr are in good agreement with 
the Lynx emission line spectrum, except for the He +2 lines 
which are predicted too weak (Table 5). Between 3.0 and 
4 Myr, WR stars appear and the He +2 problem is solved. 
Ages of 4 Myr and older produce strong discrepancies with 
the measured spectrum. We show 4 and 5 Myr models in 
Table 5 for illustration. 

Zstars = Z models 

Models with ages ^5 Myr can explain the Lynx emission 
line spectrum remarkably well, except the He +2 lines, which 
are predicted too weak. The appearance of WR stars at 2.8 
Myr solves the He +2 problem and produces a good fit to the 
observed spectrum. Ages of 6 Myr and older produce strong 
discrepancies with the measured spectrum. We only show a 
6 Myr model in Table 6 for illustration. 

To summarize, the best photoionization models imply: 

• The nebular abundances in the Lynx arc are ~10±3% 
Zq. This is in very good agreement with the measured 
abundances in §2.1.3 and with (FOSB03) within the errors. 
Lower (<7% Zq) and higher (>15% Z Q ) metallicities pro- 
duce strong discrepancies with the measured line ratios. 

• Nitrogen is overabundant. The models suggest N/O ~ 
3 x (N/0)q. This is in good agreement with the results 
obtained in §2.1.3. 

• Our models do not imply overabundance of Si relative 
to O, contrary to the results obtained by (FOSB03). The 
reason to claim Si overabundance by this authors is that 
their model predicts SiIII]/H/3 ratios a factor of ~100 below 
the measured value. Si +2 is the species found in the Lynx 
arc with the lowest ionization potential (IP). All other lines 



detected are emitted by ions with higher IP (N+ 2 and C+ 2 
are next, with 47.4 eV and 47.9 eV respectively). A hot 
black body (80 000 K in (FOSB03) models) produces a large 
supply of ionizing photons that efficiently remove Si +2 while 
still keeping a considerable amount of other species such as 
N +2 and C +2 . The same can be said about (Binette et al. 
(2003)) models. The ionization of the gas in our models is 
dominated by much colder stars that allow the survival of 
Si +2 in the nebula. 

• The electron temperature is T e ~16200±500 K, in good 
agreement within the errors with the value determined from 
the [OIII] 1663/5007 ratio (see §2.1.1). 

• Ages SS5 Myr produce good fits to the Lynx spectrum, 
except for the He + lines, which are predicted too faint in 
most cases. Only when Wolf Rayet stars appear, the He +2 
lines are properly fit. If Wolf Rayet stars are responsible for 
the ionization of He + (but see §3.1), this implies: 

- Zstars must be >5% Zq, since no Wolf Rayet stars 
are formed at lower metallicities, according to the models 
(but see §3.1). 

- the Lynx arc is undergoing a WR phase and the age 
of the burst is in the range 2.8 - 3.4 Myr, depending on 
the stellar metallicity. 



3 DISCUSSION 

3.1 The He+ 2 problem 

All models in Tables 3 to 6 with ages of ~5 Myr or younger 
(depending on the stellar metallicity) reproduce remarkably 
well most line ratios of the Lynx spectrum. Some of those 
models, however (those with no contribution from WR stars) 
cannot explain the He +2 nebular emission. 

HeIIA4686 nebular emission has been detected in the 
optical spectrum of numerous extragalactic HII regions and 
star forming galaxies (e.g. (Garnett et al. 1991)). If the 
excitation of this line is due to photoionization by hot stars, 
these must have T e ff <;55 000 K. Based on photoionization 
modeling, (FOSB03) concluded that T ef f ~80 000 K in the 
Lynx arc, in the range of values typical of Wolf Rayet stars. 

We have shown that Wolf Rayet stars are good candi- 
dates to explain the He +2 emission in the Lynx arc (they 
could also be responsible for the outflowing wind discovered 
by (FOSB03)). The non detection of the typical Wolf Rayet 
bumps found in Wolf Rayet galaxies (e.g. (Vacca & Conti 
1992)) is simply due to the fact that the stellar population 
is not detected at all (see §3.3). 

We have used two of the SEDs that best reproduce the 
Lynx spectrum (Z stars =20% Zq and 3.4 Myr; Z stars =40% 
Zq and 3.0 Myr, see §2.2) to predict the number of WR 
stars in the Lynx arc. CMK04 models predict 2.9xl0~ 4 
and 2.5 x 10~ 4 Wolf Rayet stars per solar mass for these two 
SEDs. Using the total stellar burst masses estimated from 
the H/3 luminosity (see §3.3) , we obtain 2.7xl0 4 or 2.1 xlO 4 
WR stars for the 20% Z Q and 40% Z Q SEDs respectively. 
This is in the range of values found for low redshift WR 
galaxies, where the presence of 100 to 100 000 WR stars has 
been inferred ((Vacca & Conti 1992)). 

WR stars have been proposed to explain the He + emis- 
sion in stellar ionized nebulae. However, this is an unsolved 
problem in nearby star forming objects. For this reason, we 
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Table 3. Z s tars= 0.001 = 5% Zq. Line ratios are always given relative to H/3. Discrepant values are shown in bold. The NV, NIV], 
CIV, CIII], OIIIJ1665 and [Oil] values include the fluxes of the two doublet lines. The agreement between the models with ages of 4 
Myr or less and the observations is remarkably good, except for the He+ 2 lines, which are always predicted too faint. Considering N/O 
~ (N/O)0, the models underestimate the strength of the nitrogen emission lines (see model 2 Myr a ). The problem disappears when 
nitrogen is assumed to be overabundant relative to oxygen by a factor of ~3. Models of 5 Myr and older are inconsistent with the Lynx 
arc spectrum. We only show a 5 Myr and a 6 Myr model for illustration. 
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Table 4. Z 3ta rs= 0.004 = 20% Zq. All models assume N/O ~ 3 X (N/O©). In spite of the He+ 2 problem (see text) the agreement 
between the models (age 3 Myr or less) and the observations is remarkably good. The He+ 2 problem is solved at ~3.4 Myr, when Wolf 
Rayct (WR) stars appear and ionize He+. Models of 4 Myr and older are inconsistent with the Lynx arc spectrum. We only show a 5 
Myr old model for illustration. 
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Table 5. Z s tars= 0.008 = 40% Zq. All models assume N/O ~ 3 X (N/Oq). In spite of the He+ 2 problem (sec text) the agreement 
between the models and the observations is remarkably good for an age of 3.3 Myr or less. The best models arc those with age ~3 Myr. 
when WR stars appear and ionize He+. Models of 4 Myr or more produce strong discrepancies with the observations. 
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Table 6. Z s t ars = 0.020 = Zq. All models assume N/O ~ 3 X (N/Oq). In spite of the He +2 problem (see text) the agreement between 
the models and the observations is remarkably good for ages of 5 Myr or less. The best model corresponds to 2.8 Myr of age, when Wolf 
Rayet are responsible for the ionization of He+. Models of 6 Myr or more produce strong discrepancies with the observations. 
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do not reject the possibility that another unknown ioniza- 
tion source is present. 

As an example, evidence for the existence of Wolf Rayet 
stars in the metal poor galaxy I Zw 18 (Z=l/50 Zq) has been 
reported by several authors ((Izotov et al. 1997), (Legrand 
et al. 1997)). (Stasihska & Schaerer (1999) ) showed that the 
He +2 nebular intensity can be reproduced using a radiation 
field consistent with the observed Wolf-Rayet spectral fea- 
tures in this object (see also (de Mello et al. 1998)). However, 
it is not clear whether WR stars can explain the spatial dis- 
tribution of the He +2 emission ((Vilchez & Iglesias-Paramo 
1998)). 

Izotov et al. (2004) have recently discovered the high 
ionization emission line [NeV]A3426 in the spectrum of the 
blue compact low metallicity galaxy TOL 1214-277. Fricke et 
al. (2001) detected the [FeV]A4227 line. Izotov et al. (2004) 
conclude that the stellar radiation is too soft to explain the 
existence of these high ionization species and propose that 
fast shocks or high mass X-ray binary systems could be re- 
sponsible. The possibility that ultra-luminous X-ray sources 
(ULX) are responsible for the He +2 emission in the dwarf ir- 
regular galaxy Holmberg II has been proposed by (Pakull & 
Mirioni (2002)) and (Kaaret, Ward & Zezas (2004)). ULXs 
may be be accreting normal mass (<20 Mq) black holes or 
neutron stars. An alternative possibility is that ULXs are 
high-mass X-ray binaires with super-Eddington mass trans- 
fer rates ((King et al. 2001)). 

Peimbert, Sarmiento & Fierro (1991) showed that the 
emission line spectra of giant HII regions can be altered by 
the presence of shock waves produced by stellar winds or su- 
pernova events. This is also the case in star forming galaxies 
depending on the evolutionary phase of the stellar cluster 
((Viegas, Contini & Contini 1999)). The outflowing wind 
in the Lynx arc discovered by (FOSB03) provides evidence 
for the possible presence of shocks in the Lynx arc. Such kind 
of effects are not included in our synthesis models (where the 
supernova X-ray emission is the time-average contribution 
of the events). 

A final possibility is the presence of super-soft X-ray 
sources ((Rappaport, di Stefano & Smith 1994)). One of 
the physical models proposed to explain the nature of these 
sources involves mass transfer from a main-sequence or sub- 
giant donor star to the surface of a white dwarf in a binary 
system ((van den Heuvel et al. 1992)). If this is the case in 
the Lynx arc, an additional stellar population with an age 
older than 50 Myrs must be present. 

Therefore, we do not discard the possibility that alter- 
native mechanisms to WR stars are responsible for the ion- 
ization of He + in the Lynx arc. 

3.2 The nitrogen overabundance 

We have found that the strength of the NIII] and NIV] 
lines cannot be explained assuming solar abundances. As ex- 
plained above, a possible explanation for this discrepancy is 
the assumption of constant density, rather than taking into 
account the presence of a gradient (see §2.1.3) (Luridiana et 
al. 2004, in preparation). 

An alternative explanation is that N is overabundant 
by a factor ~2-3 relative to the solar abundance (see §2.1.3 
and §2.2). Significant nitrogen excesses, although rare, have 
been found in some star forming objects. Particularly in- 



teresting are the five extremely metal poor BCGs with large 
nitrogen excesses (see (Pustilnik et al. (2004)) for a more de- 
tailed discussion on this issue). The authors propose that the 
nitrogen excesses could be a consequence of merger events 
and a short powerful starburst phase when many WR stars 
contribute to a fast enrichment of the ISM. 

The morphology ((FOSB03)) of the Lynx arc is remi- 
niscent of a merging system (several clumps A, B, C and 
D joined by a very faint arc), but this is misleading since it 
is strongly distorted by the effects of the lensing cluster. In 
their strong lensing model, (FOSB03) assume that the Lynx 
arc consists of two different clumps (A and C) of the same 
source, while B and D are mirror images of A and C re- 
spectively. A and C (or B and D) might be two components 
of a merging system. The evidence, however, is not strong 
enough to confirm this. 

3.3 The continuum problem 

We investigate here whether the starburst luminosity re- 
quired to excite the strong emission lines in the Lynx arc is 
consistent with the non detection of the stellar component. 

(FOSB03) found a discrepancy (a factor of ~20) be- 
tween the mass of the starburst needed to power the H/3 lu- 
mininosity of the Lynx arc and the mass of a starburst whose 
continuum would be just marginally detected, in agreement 
with the fact that the observed continuum is mostly nebular. 
We have done a similar exercise using the SED models that 
best reproduce the spectroscopic properties of the Lynx arc. 

Using the SEDs by (CMK04 ) we have calculated the 
stellar mass of a burst that can generate the observed Lynx 
H/3 luminosity L(Hf3)=4.3 x 10 42 erg s _1 (corrected for 
magnification, (FOSB03)). The results are shown in Table 
7 for two of the SEDs that best reproduce the emission line 
spectroscopic properties of the Lynx arc (see §2.2). 

The calculations have been done assuming covering fac- 
tor CC=1 ((FOSB03)). We obtain M siars ~9 x 10 7 M©. 
This is a lower limit, since we have assumed that all ionizing 
photons are absorbed by the gas and there is no dust red- 
dening. It is important to note that the uncertainties on the 
assumed magnification factor /i are large. (FOSB03) found 
a rather wide range of models that successfully describe the 
Lynx arc system. This introduces serious uncertainties in 
the estimated cluster mass. 

We have then calculated the stellar continuum at ~1600 
A (rest frame) emitted by these instantaneous bursts and 
compared it with the continuum level measured for the 
Lynx arc which is ~0.15 /iJy (observed frame value cor- 
rected for magnification, (FOSB03)). This corresponds to 
F ! ,=7.89xl0 -32 ergs -1 cm~ 2 Hz -1 (rest frame). The results 
are shown in Table 7. The theoretical values are a factor of 
~10 higher than the value measured for the Lynx arc. The 
discrepancy is even higher if we take the nebular continuum 
into account, which according to (FOSB03) is the dominant 
continuum component. 

As (FOSB03), we therefore conclude that the stellar 
burst that reproduces the measured H/3 luminosity would 
produce an observable continuum much (~10 times) brighter 
than that detected from the Lynx arc. Adding the nebular 
contribution, the continuum should be even brighter. The 
conclusion is not affected by uncertainties on the magnifi- 
cation factor n if the continuum and the H/3 flux scale in 
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Metallicity 


Age 


i(H/9) 




Mass n p 


Fi600 




Zstars/ Zq 


Myr 


erg s~ 


'Me' 


M Q 


erg s _1 


cm- 2 Hz" 1 


20% 


3.4 


3.2 X 


10 34 


9.4 x 10 7 


6.2x10" 


-31 


40% 


3.0 


3.5 x 


10 34 


8.6 x 10 7 


6.4x10" 


-31 


Lynx 










7.9x10" 


-32 



Table 7. Masses (Mass^) of instantaneous burst models required to reproduce the H/3 luminosity of the Lynx arc and expected 
continuum flux at 1600 A (rest frame) Fi600- The two SEDs are those that best reproduce the emission line spectrum of the Lynx arc 
(see §2.2). i(H/3) is the expected H/3 luminosity per solar mass. The continuum level expected for an instantaneous starburst of mass 
Mass 11 $ and SEDs under consideration is ~10 times higher than that measured for the Lynx arc. 



the same way with /i. A similar problem was reported by 
(Steidel et al. (2000)) regarding the giant Lya nebulae as- 
sociated with Lyman break galaxies. The authors found that 
the measured UV continuum is not enough to explain the 
high Lya luminosities of the nebulae. 

We now discuss some possible reasons for this discrep- 
ancy in the Lynx arc. 

The AGN scenario 

The previous discussion suggests that the ionizing 
source in the Lynx arc is hidden from the observers view. 
Such a geometry is characteristic of type II active galax- 
ies (AGN). (Chapman et al. (2004)) have recently suggested 
that the Lyman "blobs" found by (Steidel et al. (2000)) 
might be powered by a hidden AGN. (Binette et al. (2003)) 
suggested the possibility that the Lynx arc is a type II ac- 
tive galaxy. We, however, strongly support that this is a star 
forming object for several reasons: 

• an unusual (filtered) ionizing continuum shape is re- 
quired to explain the line ratios of the arc (Binette et al. 
2003) within the AGN scenario. Our new models using nor- 
mal (i.e. non primordial) stars show an excellent agree- 
ment with the Lynx arc measured line ratios, better than 
the filtered AGN models (see for instance, [OIII]A5007/Hb, 
NV/Hb, CIV/H/3, CIII]/H/3 in Table 1 of (Binette et al. 
(2003))). I.e. there is no need for an exotic AGN model to 
reproduce the line ratios. 

• the very narrow lines of the Lynx arc (FWHM of non 
resonant lines <100 km s _1 , (FOSB03)) are typical of star 
forming objects. AGNs usually show line widths of few hun- 
dred km s _1 . Although there are examples of extended emis- 
sion line regions in active galaxies with narrower lines, these 
cases are exceptional. 

• The strongest discrepancies of the filtered AGN model 
disappear with our new models: these are the severe un- 
derprediction of the Silll] lines and the too high electron 
temperature (~20 000 K) predicted by the filtered AGN 
models. 

We therefore reject the possibility that the Lynx arc is 
an active galaxy. 

A hidden starburst 

Hidden starbursts have been proposed to explain the 
spectroscopic properties of the metal poor blue compact 
galaxy HS 0837+4717 (Pustilnik et al. 2004). In this ob- 
ject, the extinction of the narrow emission lines from the 
giant HII region is low, but the large Balmer decrement 
(~19±4) of the broad components suggests that part of 



the current starburst is highly obscured by dust. A highly 
obscured starburst has also been revealed by mid-IR ob- 
servations of the very low metallicity galaxy SBS 0335-052 
(Z=l/41 Z Q , (Hunt, Vanzi & Thuan 2001)). Roughly 3/4 
of the star formation in this object occurs within the ob- 
scured cluster. 

It is not clear however, that this scenario solves the 
the Lynx arc problem. In HS 0837+4717 and SBS 0335-052 
the non-obscured regions emit both strong continuum and 
emission lines. I.e., if optical emission lines are detected, the 
ionizing stars are also detected. This is not the case for the 
Lynx arc. 

Differential gravitational lensing amplification 

The different spatial distribution between ionized gas 
and young stars can be the solution to the Lynx continuum 
problem. 

The discrepancy between reddening values derived us- 
ing the UV continuum and the ratios of the Balmer lines 
found in some star forming galaxies (e.g. (Mas-Hesse & 
Kunth 1998)) has been explained as due to ionized gas and 
stars having different spatial distributions and suffering dif- 
ferent extinction by dust. When this discrepancy is found, 
the UV continuum seems to be less affected by extinction 
than the emission lines. The opposite case would apply to 
the Lynx arc. Spatial decoupling between ionized gas and 
young stars has been observed in some star forming galaxies 
(e.g. (Maiz-Apellaniz et al. 1998)) and proposed by several 
authors to explain different issues (e.g. (Mas-Hesse & Kunth 
1998), (Gonzalez Delgado & Perez 2000)) 

The degree of amplification suffered by the source un- 
dergoing lensing depends strongly on both the position close 
to the critical lines and the geometry of the lensed object, so 
that separate regions can suffer different degrees of ampli- 
fication. If stars and gas have different spatial distributions 
in the Lynx arc, emission lines and stellar continuum can 
suffer different magnification due to the intervening cluster. 
Differential amplification between lines and continuum has 
been observed in several quasars (e.g. (Racine 1992)). De- 
tailed gravitational lensing models would be necessary to 
verify this possibility. 

A primordial stellar population 

The SED of metal-free stars is characterized by effective 
temperatures on the Main Sequence around 10 5 K, hotter 
than their counterparts of equal mass but finite metallic- 
ity (e.g. (Tumlinson & Shull 2000)). The UV recombina- 
tion emission lines, such as HeIIA1640, would have extreme 
equivalent widths more than an order of magnitude larger 
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than the expectation for a normal cluster of hot stars with 
the same total mass and a Salpeter IMF. Bromm, Kudritzki 
& Loeb (2001) estimated that a cluster of 10 6 Mq of popIII 
stars with a heavy IMF would produce ~2.3xl0 42 erg s _1 
in the HeIIA1640 line and a spectral luminosity per unit 
frequency of 1.8xl0 27 erg s -1 Hz" 1 . For the observed lu- 
minosity of the HeIIA1640 in the Lynx arc (~4.4xl0 42 erg 
s _1 ) a cluster of ~2x 10 6 M of primordial stars would be re- 
quired. Such a cluster would produce a spectral flux per unit 
frequency of 3.1xl0 -32 erg s" 1 cm~ 2 Hz -1 at the distance 
of the Lynx arc. This is a factor of ~2.3 below the contin- 
uum measured for the Lynx arc. The discrepancy, therefore, 
disappears and this result is consistent with the conclusion 
by (FOSB03) that the continuum is mostly nebular. 

The possibility that the Lynx arc is a popIII object was 
proposed and discussed in detail by (FOSB03). It is difficult 
to explain, however, nebular abundances as high as ~10% if 
the stellar population is primordial. 



4 SUMMARY AND CONCLUSIONS 

We have characterized the physical properties (electron tem- 
perature, density, chemical abundances) of the ionized gas 
and the ionizing stellar population in the Lynx arc, a gravi- 
tationally amplified HII galaxy at 2=3.36. The temperature 
sensitive ratio [OIII] AA1661,1666/A5007 implies an electron 
temperature T e =17300t^ K, in good agreement within the 
errors with photoionization model predictions. The UV dou- 
blets imply the existence of a density gradient in this ob- 
ject, with a highly ionized high density region (0.1-1.0 x 10 5 
cm -3 ) and a low density region (<3200 cm~ 3 ) with lower 
ionization state. 

Both the photoionization modeling and standard tech- 
niques of chemical abundance determination imply that the 
gas metallicity is ~10±3% Zq. Both methods suggest that 
nitrogen is overabundant with N/O~2.0-3x [N/0]q, unless 
a density gradient produces this apparent effect. We do not 
find evidence for Si overabundance as Fosbury et al. (2003) . 
The reason is the different shape of the ionizing continuum 
assumed by these authors (80 000 K black body, much hotter 
than the dominant ionizing stars in our models). 

Photoionization models imply that the ionizing stars 
have very young ages ^5 Myr. Since the emission lines trace 
the properties of the present burst only, nothing can be said 
about the possible presence of an underlying old stellar pop- 
ulation. Instantaneous burst models with Z siar >5% Z and 
ages ~2.8-3.4 Myr (depending on Z s t ar ), are in excellent 
agreement with the Lynx spectrum, including the strong 
He+ 2 emission. At this age Wolf Rayet stars make an impor- 
tant contribution to the hard ionizing luminosity and they 
are responsible for the excitation of the He +2 emission. In 
such case, we infer the existence of ~2.5xl0 4 WR stars in 
the Lynx arc. Alternative excitation mechanisms for He +2 , 
however, cannot be discarded. 

Therefore, the Lynx arc is a low metallicity HII galaxy 
that is undergoing a burst of star formation of ^5 Myr age. 
One possible scenario that explains the emission line spec- 
trum of the Lynx arc, the strength of the nitrogen lines and 
the strong He +2 emission is that the object has experienced 
a merger event that has triggered a powerful starburst phase. 



Wolf Rayet stars have been formed and contribute to a fast 
chemical enrichment of the interstellar medium. 

As Fosbury et al. (2003), we find a factor of >10 discrep- 
ancy between the mass of the instantaneous burst implied 
by the luminosity of the H/3 line and the mass implied by 
the continuum level measured for the Lynx arc. We have dis- 
cussed several possible solutions to this problem. The most 
satisfactory explanation is that gas and stars have different 
spatial distribution so that the emission lines and the stellar 
continuum suffer different gravitational amplification by the 
intervening cluster. Detailed gravitional lensing models are 
needed to test the vailidity of this scenario. 
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